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ABSTRACT 


This  thesis  provides  a  method  to  utilize  existing  digital 
simulations  and  missile  reliability  data  to  determine  the 
probability  of  kill  (Pk)  versus  range  for  the  Seasparrow  for  each 
system  state.  A  procedure  is  presented  that  uses  the  Pk  versus 
range  curve  and  the  firing  doctrine  to  detearmine  the  probability  of 
defeating  ein  antiship  cruise  missile  attack  for  input  into  a 
baseline  system  effectiveness  model.  A  method  to  optimize  the 
firing  doctrine  using  linear  prograutiming  techniques  is  presented. 
Implemented  using  the  General  Algebraic  Modeling  System  (GAMS)  ,  the 
optimized  firing  doctrine  provides  a  significant  enhancement  to  the 
probability  of  defeating  an  ASCM  attack.  The  thesis  enhances  the 
Baseline  Systems  Effectiveness  model  for  the  NATO  Seasparrow 
Missile  System  (NSSMS)  developed  by  Professor  W.M.  Woods  of  the 
Naval  Postgraduate  School  for  the  Naval  Warfare  Assessment  Center, 
Corona,  California. 
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The  reader  is  cautioned  that  con^uter  programs  developed  in 
this  research  may  not  have  been  exercised  for  all  cases  of 
interest.  While  every  effort  has  been  made,  within  the  time 
available,  to  ensure  that  programs  are  free  of  con^utational  and 
logic  errors,  they  cannot  be  considered  validated.  Any  application 
of  these  programs  without  additional  verification  is  at  the  risk  of 
the  user. 
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I  ZBITRODOCTION 


The  NATO  Seasparrow  Missile  System  (NSSMS)  is  the  primary 
point  defense  system  for  the  United  States  Navy.  Like  laany 
modem  weapon  systems,  the  Seasparrow  is  designed  with  a 
variety  of  alternate  operating  modes  for  use  in  the  event  of 
coR^onent  failures.  Although  this  is  desirable  from  an 
operational  standpoint,  it  creates  difficulties  in  assessing 
the  utility  of  the  system.  The  Navy  has  traditionally  tracked 
system  availability  by  recording  what  percentage  of  the  time 
a  system  is  fully  operational,  usually  referred  to  as  A„.  This 
method  has  several  drawbacks  when  applied  to  the  NSSMS.  The 
major  drawback  is  the  lack  of  accounting  for  degraded  modes  of 
operations.  Additionally,  it  provides  no  measure  of  the 
sy*- terns  capability,  only  reliability  and  maintainability.  To 
address  these  concerns,  system  effectiveness  models  were 
developed.  In  essence,  a  system  effectiveness  model  svons,  for 
all  operating  states,  the  product  of  the  expected  proportion 
of  time  a  system  will  be  in  a  given  state  (system  state 
availability) ,  the  probability  the  system  transitions  to  a 
different  state  during  an  engagement,  and  the  capability  of 
the  system  for  the  system  state  in  which  the  system  is 
operating  at  the  time  of  missile  launch.  This  paper  provides 
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an  analytical  framework  to  determine  the  probability  of  the 
NSSMS  defeating  an  antiship  cruise  missile  (ASCM)  attack,  for 
input  into  a  NSSMS  system  effectiveness  model  developed  for 
the  Naval  Warfare  Assessment  Center,  Corona,  CA  (NWAC) .  This 
probability  is  one  factor  in  the  NSSMS  system  effectiveness 
equations . 

A.  PURPOSE 

A  method  is  developed  to  determine  the  probability  of  the 
NSSMS  defeating  an  antiship  cruise  missile  (ASCM)  attack  based 
on  the  operating  state  of  the  NSSMS  at  the  time  the  missiles 
are  launched.  The  thesis  will  explore  the  use  of  existing 
digit?'!  simulations  developed  by  the  Naval  Air  Warfare  Center, 
China  Lake,  CA,  (NAWC)  to  develop  an  estimate  of  the  single 
shot  probability  of  kill  (Pk)  for  the  Seasparrow  when  the 
NSSMS  is  in  system  state  S^  at  the  time  of  launch.  This 
estimate  will  be  combined  with  missile  reliability  data  and 
used  in  conjunction  with  the  firing  doctrine  to  evaluate  the 
overall  probedsility  of  kill  (Pk)  for  the  NSSMS  against  a 
representative  subsonic  ASCM.  Finally,  a  method  to  use  linear 
progreunming  techniques  to  optimize  the  firing  doctrine  versus 
a  single  target  will  be  exeunined. 


2 


II 


BACKOROUND 


In  an  effort  to  improve  the  way  in  which  the  reliability 
of  the  NATO  Seasparrow  Missile  System  (NSSMS)  is  tracked,  the 
Naval  Warfare  Assessment  Center  (NWAC)  funded  Professor  W.M. 
Woods  of  the  Naval  Postgraduate  School  to  develop  a  baseline 
systems  effectiveness  model  for  the  system  [Ref.  1]  .  The  model 
he  developed  is  based  on  methodology  developed  by  the  Weapon 
Systems  Effectiveness  Industrial  Advisory  Committee  [Ref.  2]. 
The  baseline  model  recognizes  seven  system  states,  S^  to  S7, 
and  is  written  as 

7  7 

^  (2.1) 

where  AS^  is  the  operational  availability  of  system  state  S^, 
DS^j  is  the  probability  the  system  transitions  from  state  S^ 
to  state  during  an  engagement,  and  CSj  is  the  capability  of 
the  system  when  in  state  S j .  This  paper  provides  a  method  to 
determine  the  value  of  CSj  using  digital  simulation  and 
missile  reliability  data. 

CSj  is  defined  as  the  capability  of  the  system  to 
defeat  the  chosen  threat  given  eight  missiles  are  loaded  at 
the  beginning  of  the  scenario  [Ref.  1] .  In  the  baseline  model, 
CSj  is  defined  as  the  probability  of  success  given  that  two 
missiles  are  fired  at  the  inbound  target.  This  is  written  as 
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CSj=Pjl)Pkj{l)  *Pj2)Pkj(2) 


(2.2) 


where  PlIh)  is  the  probability  of  successfully  firing  n 
missiles  and  Pk^  (n)  is  the  probability  of  successfully 
destr'jying  the  target  given  that  n  missiles  were  fired  in 
system  state  Ss .  It  should  be  noted  that  equation  (2.2) 
incorporates  the  missile  reliability  term,  R^,,  in  the  original 
CSj  formula  into  the  Pkj  (n)  term  for  brevity.  The  formulation 
and  calculation  of  CSj  will  be  further  discussed  in  Chapters 
IV  and  V. 

A.  SYSTEM  CC»lFIGimATZC»I 

The  NSSMS  comes  in  two  basic  configurations:  one  fire 
control  direc..or  or  two  fire  control  directors.  Each  director 
has  its  own  transmitter,  receiver,  signal  data  processor 
(SDP) .  and  radar  system  controller  (RSC) .  An  eight  cell  box 
launcher  is  used  in  both  configurations  for  launching  the 
missiles.  For  purposes  of  systems  effectiveness,  if  a  two 
director  configured  ship  has  a  casualty  in  only  one  director, 
the  system  is  considered  to  be  fully  functional. 

B.  SYSTEM  STATES 

The  baseline  system  effectiveness  model  for  the  NSSMS 
recognizes  seven  different  system  states.  Although  there  are 
more  than  seven  possible  degradations  to  the  system,  only 
seven  were  considered  to  occur  with  sufficient  frequency,  or 
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have  a  sufficient  impact  on  probability  of  success,  to  justify 
inclusion  in  the  model.  The  system  states  and  their  associated 
equipment  casualties  are  listed  in  Table  1  [Ref.  l,p.  6]: 

Table  It  NSSMS  OPKRATIHQ  STATES 


System 

State 

Operating 

Mode 

Equipment  Casualty 

Si 

Normal 

All  up 

S2 

Normal 

All  up  except  fault  code  100: 

Transmitter  power  output  degraded. 

S3 

Normal 

All  up  except  fault  code  106:  Signal 

Data  Processor  range  signal  degraded. 

S4 

Normal 

All  up  except  fault  code  112:  Range 
deviation  too  large. 

S5 

Normal 

All  up  except  POC/RSC  which  has  only 
co'trol  power  sufficient  to  assign, 
tune,  and  fire  missile  in  a  normal 
manner. 

S6 

Cl  (Slaved 
Illuminator) 
or  LLLTV 

All  up  except  SDP  or  receiver. 

S7 

Computer 

complex 

Casualty 

(CCC) 

All  up  except  digital  computer  or 
signal  data  converter  (SDC) . 

C.  MZSSZUE  COMMANDS 

The  RIM-7  Seasparrow  missile  is  a  derivative  of  the  AIM-7 
Sparrow  Air-to-Air  missile.  It  has  been  optimized  for  the 
point  defense  role  for  use  against  antiship  cruise  missiles. 
The  missile  employs  semi-active  radar  homing  zmd  pre-launch 
missile  commands  cannot  be  updated  in  flight.  Degraded  missile 
commands  are  the  primary  effect  of  equipment  degradations  in 
system  states  83,84,85,  and  8^.  The  seven  missile  commands  are 
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listed  below  with  one  system  command,  launcher 
super-elevation,  which  is  also  affected  by  system  degradations 
that  may  affect  missile  performance. 

a.  Seeker  Antenna  (Head:  Aim (pitch  and  yaw) 

b.  English  Bias 

c.  Sweep  Select 

d.  Sweep  Control 

e.  Video  Inject 

f .  Range  at  Launch 

g.  Channel  Designate/HOH 

h.  Launcher  Super-Elevation 

Reference  3  provides  a  complete  discussion  of  these  missile 
commands.  In  the  event  of  a  system  casualty  that  prevents  the 
fire  control  system  from  providing  inputs  to  the  missile, 
state  Sg  for  exeui^le,  default  values  are  utilized.  Of  these 
commands.  Range  at  Launch  is  evaluated  as  having  the  greatest 
effect  on  missile  performance. 

D.  DIGITAL  SIMDLATICmS 

Due  to  the  limited  number  of  actual  firings,  it  is 
necessary  to  conduct  simulations  to  determine  the  missiles 
probability  of  success.  The  NATO  Seasparrow  program  office 
(NSPO)  ,  maintains  a  network  of  simulations  which  are  available 
throughout  the  Seasparrow  technical  community  [Ref.  4].  The 
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Navy  organization  responsible  for  the  Seasparrow  simulations 
is  the  Naval  Air  Warfare  Center,  China  Lake,  CA  (NAWC) .  NAWC 
maintains  two  Seasparrow  flight  simulations  that  offer 
sufficiently  accurate  modelling  to  give  valid  results:  a  six 
degree  of  freedom  digital  simulation  and  a  six  degree  of 
freedom,  hardware-in-the-loop  simulation.  Due  to  the  order  of 
magnitude  different  in  cost  between  the  digital  simulation  and 
the  hardware-in-the-loop  simulation,  the  digital  simulation  is 
the  preferred  choice.  In  addition  to  the  flight  simulation, 
NAWC  maintains  several  endgeume  lethality  simulations  that  will 
be  utilized,  including  a  target  detecting  device  (TDD) 
simulation,  a  warhead  fragmentation  pattern  simulation,  and  a 
target  vulnerability  model. 

The  four  simulations  are  used  sequentially  to  determine 
the  probability  of  success  for  each  simulation  trial.  The 
flight  simulation  is  used  to  determine  the  flight  path  of  the 
missile  from  launch  to  intercept.  The  last  few  seconds  of  the 
missile  and  target  flight  paths  are  then  entered  into  the  TDD 
simulation  to  calculate  the  warhead  detonation  point.  After 
determining  the  warhead  detonation  point,  the  relative 
positions  and  velocities  of  the  missile  and  target  are  fed 
into  the  warhead  simulation  to  determine  how  many  fragments 
struck  the  target.  The  estimated  fragmentation  damage  to  the 
target  is  then  compared  with  the  target  vulnerability  model  to 
determine  the  probability  of  the  target  being  destroyed. 
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III.  SIMDIATIOM  PLAN 

Due  to  the  limited  number  of  actual  missile  firings  versus 
ASCHs,  simulations  are  required  to  determine  the  intact  of 
range  and  system  degradations  on  the  probability  of  kill  (Pk) 
for  the  Seasparrow  in  each  system  state.  Digital  simulations 
offer  the  advantage  of  s^abstantially  lower  cost  and 
significantly  shorter  turn  around  time,  when  compared  to 
hardware- in- the- loop  simulations. 

In  any  systems  effectiveness  model,  a  specific  scenario  is 
required  to  determine  the  probability  of  success  of  the 
system.  The  baseline  systems  effectiveness  model  [Ref.  1] 
provides  three  Seasparrow  missiles  configurations  and  three 
possible  ASCMs  to  be  evaluated.  The  simulation  plan  in  this 
thesis  is  developed  for  the  RIM-7P  configuration  of  the 
Seasparrow  versus  the  Ebcocet  cruise  missile.  This  selection 
was  chosen  since  all  elements  required  (flight  simulation,  TDD 
simulation,  warhead  simulation,  and  target  vulnerability 
model)  are  readily  available  at  NAWC  China  Lake  [Ref.  5]. 

A.  BACKGROUND 

Initially,  it  was  assumed  that  the  probability  of  kill  for 
the  Seasparrow  could  be  represented  by  a  single  value  for  the 
entire  range  of  the  missile.  Conversations  with  the  technical 
coiMCiunity  [Refs.  5&6]  indicated  that  this  would  be  an  over 
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simplification.  This  would  be  especially  true  of  any  degraded 
modes  of  operation  for  the  system-  It  is  expected  that  Pk  will 
increase  as  the  velocity  of  the  Seasparrow  increases,  and 
decrease  when  the  velocity  of  the  missile  decreases.  This 
effect  should  be  especially  noticeable  at  or  near  maximum 
range  when  the  missile  is  coasting  and  has  very  little  energy 
available  for  terminal  maneuvers.  It  is  anticipated  that  the 
Pk  versus  range  curve  will  look  roughly  as  shown  in  Figure  1. 


Figure  1<  KMpected  Pk  vm  Range  curve 


It  should  be  noted  the  values  for  Pk  indicated  above  are 
for  illustrative  purposes  only  and  do  not  represent  the  actual 
probability  of  kill  for  the  Seasparrow  missile. 
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B.  RBQnZRXMBNTS 

The  simulation  plan  must  conform  to  several  constraints. 
The  plan  must  cover  the  full  range  of  the  Seasparrow,  have  a 
sufficient  number  of  test  points  to  provide  good  shape 
definition  for  the  PK  versus  Range  curve,  have  enough  trials 
at  each  point  to  provide  a  reasonable  confidence  interval  for 
Pk,  and  meet  the  cost  constraints  discussed  in  Reference  1. 

C.  ASSDMPTIOHS 

As  a  result  of  limitations  of  the  digital  simulations, 
several  assvmqptions  were  made  when  the  simulation  plan  was 
created.  First,  the  engagement  is  in  an  electronic  counter¬ 
measures  (ECM)  clear  environment,  at  least  in  the  spectrum 
covered  by  the  NSSMS  director.  The  engagement  is  limited  to  a 
single  inbound  ASCM,  and  each  missile  is  assumed  to  act 
independently  from  all  others.  Finally,  the  firing  ship  is  the 
target  for  the  ASCM,  euid  the  Exocet  flies  its  optimum  profile, 
detecting  the  firing  ship  at  maximum  range  and  commences  sea 
skimming  as  soon  as  possible. 

D.  ErySTSM  STATSS 

Due  to  the  scenario  chosen,  RIM-7P  versus  Exocet  in  a  ECM 
clear  environment,  the  seven  system  states  can  be  reduced  to 
five  states  by  combining  states  S5  and  S7.  Degradations  to 
the  FOC,  S5,  are  only  important  in  an  ECM  environment,  while 
a  con^uter  or  signal  data  converter  failure,  S7,  only  affects 
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the  ability  to  acquire  the  target  and  should  have  negligible 
effect  on  the  probability  of  success. 

E.  SZMOLATZCEI  PLAN 

To  explore  the  Pk  versus  reuige  dependency,  thirteen 
intercept  points  were  chosen,  with  the  points  spaced  closest 
at  the  extremes  of  range.  Intercepts  closer  than  five  percent 
of  the  mzucimum  range  are  not  evaluated  as  this  is  inside  the 
remge  of  the  close  in  weapons  systems.  A  semiple  size  of  fifty 
was  adopted  as  a  conpromise  between  cost  and  confidence  limit 
precision.  Complete  details  of  the  Simulation  Plan  can  be 
found  in  Appendix  A. 
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ZV.  mkTA  MIALYSZS 


Once  the  simulations  are  con^leted,  the  resultant  data  and 
the  missile's  reliability  data  needs  to  be  analyzed  and 
combined  to  determine  the  Pk  versus  range  values.  Once  these 
values  are  determined,  they  can  be  con^ared  to  the  firing 
doctrine  to  find  Pk(n)j,  the  probability  of  killing  the  target 
given  n  missiles  are  successfully  fired  in  system  state  Sj, 
and  PL(n),  the  probedjility  of  successfully  firing  n  missiles 
during  an  engagement. 

A.  KSTIMATINO  PROBABILITY  OF  KILL 

An  estimate  of  the  probability  of  kill  for  any  single 
missile  is  given  by  conibining  the  probability  that  the  missile 
guides  correctly  once  it  la\mches,  R^q,  the  probability  that 
the  missile  doesn't  fail  in  flight,  R^,  and  the  probability 
of  success  given  that  first  two  conditions  are  met,  Pkg.  The 
reliability  information  can  be  derived  from  missile  flight 
data  available  from  the  flight  analysis  department  at  NWAC, 
while  the  probability  of  success  is  derived  from  the 
simulations. 

1 .  Detezmlning  Pkg 

The  fifty  values  at  each  simulation  point  constitute 
fifty  independently,  identically  distributed  (IlD)  estimates. 


I>ksi»  i=l,2,  .  .  .  ,  50,  of  the  probability  of  kill,  Pkg.  By  the 
Central  Limit  theorem,  the  mean,  x,  of  the  fifty  estimates 
should  have  a  nearly  normal  distribution  with  mean  Pkg.  If  a 
denotes  the  standard  deviation  of  the  probability  distribution 
for  Pkg,  then  the  quantity 


(x-Pkg) 

(a/Vn) 


(4.1) 


has  a  nearly  standard  normal  distribution.  Since  the  standard 
deviation,  o,  is  not  known,  the  seutple  standard  deviation,  s, 
is  substituted  for  o  in  equation  (4.1).  The  probability 
distribution  of  the  resulting  expression  is  approximated  by  a 
t-distribution  with  n-1  degrees  of  freedom.  The  resulting 
formula  for  a  95%  lower  confidence  limit,  ^*ks(^^),  for  Pkg  is 


(4.2) 


where 


(4.3) 


(4.4) 

and  n  is  the  number  of  trials  (replications  of  simulations) 
and  t(n-i),o.95  95th  percentile  of  a  t  probability 
distribution  with  n-1  degrees  of  freedom  . 
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2.  Mlsslltt  RAllabillty 

After  Pkg  is  determined,  it  must  be  multiplied  by  the 
two  missile  reliability  terms  to  determine  the  value  for  Pk  at 
any  intercept  range.  The  first  term,  the  probability  that  the 
missile  guides  once  it  is  launched,  R^^q,  is  singly  the 
probability  that  all  components  of  the  missile  function 
properly  at  launch.  The  second  term,  the  probability  of  no 
in-flight  failures,  R^p,  varies  with  the  duration  of  the 
flight  [Ref.  7].  Missile  reliability  can  be  adequately 
represented  by  an  exponential  life  distribution  reliability 
model,  where  the  probability  of  a  successful  flight  is  defined 
as 


P {Reliable  flight)  (4.5) 

where  R[,q  is  the  probability  that  all  missile  conponents  are 
functioning  correctly  at  launch  and  R^p  is  the  probability 
that  a  missile  functions  correctly  until  intercept.  R)Qr  has  an 
exponential  life  distribution. 

a.  SmtiauLting 

A  point  estimate,  of  R^  the  probability  of  the 
missile  to  correctly  guide  once  launched,  is  the  ratio  of 
guidance  successes  divided  by  the  total  number  of  launches.  A 
guidance  failure  is  defined  as  a  missile  that  launched  but 
failed  to  guide  properly  because  not  all  coirponents  functioned 
at  launch.  The  nvunber  of  successes  is  the  number  of  tests 
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minus  the  number  of  guidance  failures.  A  100(l-a)%  lower 
confidence  limit,  ^mgcd  '  ^mg  can  be  obtained  by  looking  up 

the  appropriate  value  in  a  table  of  confidence  limits  for 
binomial  attribute  data  using  existing  tables  such  as 
Reference  8.  Alternatively, 

where  n  is  the  number  of  launches,  f  is  the  number  of 
failures,  and  is  the  100(l-o)%  percentile  point  of  the 

F-distribution  with  2(n-f-»-l)  and  2f  degrees  of  freedom.  If 
RIM-7M  and  RIM-7P  missiles  have  sufficient  commonality  of 
con^onents,  the  number  of  test  firings  could  be  pooled  for 
determining  ftfjg  and  Amgcd* 
b.  SmtiMting 

The  second  component  of  missile  reliability,  Rjjp, 
involves  failures  that  occur  in  flight.  Due  to  the  limited 
number  of  failures  that  could  be  attributed  to  in-flight 
failures,  as  opposed  to  undiscovered  failures  that  occurred 
prior  to  launch,  a  qualitative  statement  on  the  distribution 
of  in-flight  failures  is  inpossible.  For  sinplicity,  the 
failure  rate.  A,  is  assumed  to  be  constant,  which  inplies  an 
e3q>onential  life  distribution  for  the  in-flight  time  to 
failure.  The  resultant  term  for  in-flight  reliability  is 

(4.7) 

and  its  estimate  is 
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(4.8) 


where 


t  is  the  time  of  flight  to  the  intercept  point,  f^j  is  the 
number  of  in-flight  failures,  and  Tj  is  the  total  flight  time 
for  all  missile  test  firings.  Again,  if  the  RIM-7M  and  RIM-7P 
have  sufficient  commonality,  the  test  results  can  be  pooled. 

3 .  Estlaatlna 

The  estimate  of  the  probability  of  kill  derived  from 
the  simulations  is  multiplied  by  the  estimate  of  the 
probability  that  the  missile  guides  correctly,  and  the 
estimate  of  the  missiles  in-flight  reliability  for  each  of  the 
thirteen  san^le  points;  that  is 

i=l,2,...,13  (4.10) 

where  ftup  is  recalculated  for  each  intercept  point.  The  result 
is  an  estimate  for  Pk  at  each  of  the  thirteen  intercept 
ranges.  It  is  then  necessary  to  fit  a  curve  to  the  points  to 
allow  estimation  of  Pk  at  any  possible  intercept  range. 

4.  Curve  Fitting 

There  are  variety  of  commercially  available  graphics 
programs  that  perform  curve  fitting.  Given  the  expected  shape 
of  the  Pk  versus  range  curve,  a  high  order,  seventh  to  tenth. 
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polynomial  curve  should  best  fit  the  data.  Figure  2  below 
shows  a  ninth  order  polynomial  curve  fit  generated  by  the 
cor  mercially  available  program  DeltaGraph  Professional,  but 
many  other  graphics  programs  support  curve  fitting. 


Pk  VS  Intercept  Range 

NInlh  Order  Polynornal  (Xive  Fit 
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Figure  2x  Sample  Curve  Fitting 


When  the  best  fit  has  been  determined  for  the  data, 
the  formula  used  to  generate  the  curve  provides  an  estimate 
for  Pk  for  any  intercept  range. 

B.  X8TZMATZ1KI  Rn, 

Once  the  value  of  l*k  have  been  calculated  for  all  possible 
intercept  ranges,  the  next  component  required  to  calculate  CSj 
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is  PlIh),  the  probability  of  being  able  to  launch  n  missiles. 
The  factor  required  to  calculate  P^in)  is  Rjjl,  the  probability 
that  the  missile  successfully  launches.  A  success  is  defined 
as  a  missile  that,  after  the  firing  signal  was  generated, 
launched.  A  failure  is  when  the  launch  signal  is  generated, 
and  the  missile  does  not  leave  the  launch  rail.  Since  the 
NSSMS  detects  misfires  and  automatically  atteirpts  to  fire  the 
next  available  missile,  these  failures  should  not 
significantly  affect  the  probability  of  kill,  but  do  impact 
the  probability  of  being  able  to  launch  a  specified  number  (  : 
missiles . 

A  point  estimate,  of  the  probability  of  the 
missile  to  correctly  launch,  is  the  ratio  of  the  number  of 
successful  launches  divided  by  the  total  number  of  attempted 
la\anches.  A  100(l-a)%  lower  confidence  limit,  ftjiLd,) ,  for  ^ 
can  be  obtained  by  using  tables,  as  for  ftjjccD  •  Alternatively, 


6 

nor  +  (j3AL--nsr+l)F,. 


(4.11) 


where  ngL  is  the  number  of  successful  launches,  n^  is  the 
number  of  attempted  launches,  and  is  the  100(1-0)% 

percentile  point  of  the  F-distribution  with  2(nsL+l)  and 
2  (n^u^-ngL)  degrees  of  freedom.  As  with  the  other  reliability 
terms,  if  the  RIM-7M  and  RIM-7P  missiles  have  sufficient 
commonality  of  components,  the  number  of  attenpted  firings 
could  be  pooled  for  determining  and  • 
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V.  CAPABILITY  OF  SYSTEM 


Once  the  estimates  for  the  value  for  Pk  at  any  intercept 
remge,  and  the  relieibility  of  the  missile  to  launch, 

have  been  made,  it  is  possible  to  estimate  CSj,  the  capability 
of  the  system  in  state  Sj.  In  the  baseline  model  CSj  was 
defined  as  the  probedsility  of  success  given  two  missiles  were 
fired.  Since  the  NSSMS  has  sufficient  range  to  allow  a  shoot- 
shoot  -look-  shoot  -shoot  doctrine  against  most  threats,  CSj  has 
been  redefined  as  the  probability  of  success  given  up  to  four 
missiles  are  fired  during  an  engagement.  This  change  allows 
CSj  to  more  accurately  reflect  the  expected  en^loyment  of  the 
system  versus  subsonic  cruise  missiles.  CSj  is  made  up  of  two 
conponents,  Pk(n)j,  the  probability  of  success  given  that  n 
missiles  are  fired  in  state  Sj,  and  Pi^(n),  the  probability  of 
successfully  firing  n  missiles  given  eight  missiles  are 
available  at  the  beginning  of  the  engagement.  The  new  equation 
for  CSj  is  written  as 


CS^=Y^P^{n)Pk(n) 

ll*l 


(5.1) 


where  n  is  the  number  of  missiles  fired  duxxug  the  engagement. 
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A.  ISTZNATIMO  Pi.(n) 

To  estimate  the  probability  of  firing  one,  two,  three,  or 
four  missiles,  the  missile  launch  reliability  term,  is 

required.  For  the  case  of  one,  two,  or  three  missiles  fired, 
the  binoitdal  formula  provides  the  probability  of  each 
occurrence . 

JIS3.  (5.2) 

^l(4)  includes  the  probability  that  four  or  more  missiles  are 
edsle  to  be  fired,  and  is  written  as 


®  (5.3) 

To  sia^lify  the  calculations,  f*L(4)  can  be  thought  of  as  one 
minus  the  probability  that  n  is  less  than  four,  which  is 
written  as 

i>i(4)=l-((2^i:(n))+(l-^«)*)  ■  (5-4) 

o«i 

Once  each  term  for  Pli{n)  has  been  estimated,  it  is  multiplied 
by  the  appropriate  Pk(n)j  term  to  find  CSj . 
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B.  BSTZMXTIliO  Pk(n)j 

To  estimate  the  value  of  Pk(n)j,  the  firing  doctrine  is 
reviewed  to  determine  at  what  range  the  intercepts  would 
occur,  and  the  appropriate  values  of  Pki  are  determined  from 
the  Pk  versus  range  curve.  The  formula  for  Pk(n)j  is  then 

(5.5) 


with  Rli  being  the  estimated  range  at  intercept  for  missile  i. 
Once  the  value  of  f>k(n)j  has  been  determined  for  all  values  of 
n,  the  values  can  be  multiplied  by  their  respective 
terms  and  summed  to  find  CSj .  When  estimating  ^>k(n)j  each 
missile  is  assumed  to  act  independently,  and  their  is  no 
acciamulation  of  damage  to  the  target.  This  sin?)lifies  the 
calculation  and  provides  the  most  conservative  estimate  of 
success . 

For  system  state  Sg,  SI/LLLTV  mode,  on  those  ships  that  do 
not  have  a  three  dimensional  radar  to  provide  positioning  data 
to  the  NSSMS  director,  replenishment  ships  for  exeunple,  the 
low  light  level  television  is  used  for  tracking,  and  the 
maximum  range  of  detection,  and  thus  intercept,  must  be 
shortened  accordingly.  Specifications  for  LLLTV  expected 
detection  ranges  can  be  found  in  Reference  9 . 
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C.  SKJOtCH  RADARS 

The  detection  range  plays  an  iir^ortant  role  in  detenaining 
the  capability  of  the  system.  This  is  especially  true  for  the 
Coir¥>uter  Coir^lex  Casualty  mode,  state  S-j,  which  relies  on 
manual  acquisition  of  the  target  with  the  NSSM5  fire  control 
radar.  Adding  search  radar  performance  and  acquisition  times 
to  the  baseline  system  effectiveness  model  would  enhance  the 
accuracy  of  the  baseline  model,  but  would  entail  a  substantial 
increase  in  the  complexity  of  the  model.  Although  not  covered 
in  detail,  a  method  to  incorporate  search  radar  performance  is 
provided  as  a  basis  for  future  enhancements  to  the  baseline 
system  effectiveness  model. 

1.  Modlflcatloiui  to  BasollzMit  Modal 

The  baseline  system  effectiveness  model  is  modified  to 
account  for  search  radar  performance  and  reliability  by 
introducing  an  additional  factor  in  each  term  of  the  original 
formula.  The  system  effectiveness  equation  becomes 

*  7  7 

^  ^  (5.7) 

where  ACj^  is  the  availability  of  each  shipboard  search  radar 
configuration  Cy^,  m  is  the  number  of  possible  configurations 
of  the  search  radars,  and  is  the  capability  to  defeat 

the  threat  given  that  the  NSSMS  is  in  state  Sj,  and  the  search 
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radars  are  in  configuration  C]^.  The  other  terms  remain 
unchanged  from  the  baseline  model. 

Each  search  radar  configuration,  Cj^,  is  determined  by 
the  search  radar  configuration  of  the  ship.  A  Spruance  class 
destroyer  might  have  eight  configurations  as  shown  in  table  2. 
A  Wichita  class  replenishment  oiler,  AOR,  might  have  only  two 
configurations,  MK23  Target  Acquisition  System  (TAS)  up  and 
TAS  down. 


Table  2t  DD963  SEARCH  RADAR  CONFZODRATZONS 


Ck 

Radar  Status 

Cl 

TAS  up,  SPS-40  up,  SPQ-9  up 

C2 

TAS  down,  SPS-40  up,  SPQ-9  up 

C3 

TAS  up,  SPS-40  down,  SPQ-9  up 

C4 

TAS  down,  SPS-40  down,  SPQ-9  up 

C5 

TAS  up,  SPS-40  up,  SPQ-9  down 

C6 

TAS  down,  SPS-40  up,  SPQ-9  down 

C7 

TAS  up,  SPS-40  down,  SPQ-9  down 

Ce 

TAS  down,  SPS-40  down,  SPQ-9  down 

Depending  on  the  search  radar  in  question,  there  could 
be  several  possible  operating  states  that  could  effect  the 
capability  of  the  NSSMS  to  engage  the  target.  For  instance,  on 
some  ships,  the  TAS  computer  links  the  combat  direction  system 
and  the  NSSMS  together  (Refs.  9&10]  .  If  the  TAS  con^uter 
fails,  the  NSSMS  operators  must  manually  acquire  all  targets. 
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increasing  the  acquisition  time  significantly.  This  would 
result  in  three  states  for  TAS,  all  up,  radar  down  but 
computer  up,  and  all  down,  instead  of  two,  and  the  number  of 
configurations  would  increase  accordingly. 

The  new  term  for  the  capability  of  the  system,  CSj 

needs  to  incorporate  three  additional  conponents;  the  expected 
range  of  first  detection,  the  acquisition  time,  and  the 
ejqjected  probability  of  success  given  the  target  is  acquired 
at  a  specific  range.  The  addition  of  these  terms  would  allow 
the  effects  of  system  states  S^,  FOC/RSC  degraded,  and  S7,  CCC 
mode,  to  be  isolated  from  system  state  S^,  all  up,  resulting 
in  a  better  eetimate  of  the  capability  of  the  system. 
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VI.  OPTIMIZING  FIRING  DOCTRINE 


Once  the  Pk  versus  range  curve  is  known,  it  is  possible  to 
optimize  the  firing  doctrine  to  maximize  the  probability  of 
defeating  the  threat.  Ideally,  if  the  optimization  were  a 
linear  program,  one  could  take  advantage  of  commercially 
available  optimization  software  such  as  GAMS  or  LINDO.  Pk(n) 
is  maximized  by  minimizing  the  probability  of  missing.  The 
probability  of  missing  can  be  written  as 


4 


(6.1) 


Although  the  probability  of  missing  is  still  a  non-linear 
equation,  by  taking  its  natural  logarithm,  it  becomes  a  linear 
equation. 


4 


ln{l-(PJcJJ?Ji))  , 


(6.2) 


which  is  readily  solvable  as  a  modified  set  covering  problem. 

A.  LINEAR  PROGRAM 

The  first  requirement  to  allow  the  use  of  linear 
programming  techniques  is  to  break  up  the  Pk  vs  range  curve 
into  discrete  time  intervals  with  constant  Pk  within  the 
interval.  The  use  of  one  second  intervals  based  on  the  time  of 
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first  possible  launch  (maximum  range  intercept)  until  the  last 
possible  launch  (minimum  range  intercept)  results  in 
approximately  one  hundred  and  twenty  time  intervals.  This 
provides  sufficient  resolution  to  minimize  the  effects  of 
discretizing  the  continuous  data,  but  makes  the  problem  small 
enough  to  be  readily  solvable.  For  each  launch  time  interval, 
the  appropriate  intercept  range  as  the  percentage  of  loaximum 
range  is  entered  into  the  formula  selected  for  the  curve 
fitting  to  return  the  Pk  for  the  time  interval.  Solving  the 
following  formulation  results  in  the  optimum  firing  doctrine. 

1.  Zzidex 

t;  The  time  interval  in  which  the  missile  is  fired,  t  is 
indexed  from  1,  the  mzucimum  range  intercept  lavmch,  to 
^max'  the  minimum  range  intercept  launch. 

u:  An  alias  for  t,  u  offers  the  same  index  for  use  with 
matrix  notation. 

2.  Data 

P(.:  The  natural  log  of  the  probability  of  missing  (l-Pkj) 

for  a  missile  fired  during  time  interval  t. 

A  confounding  matrix,  a  one  indicates  a  missile 
launched  at  time  interval  u  will  be  airborne  during  tie 
interval  t.  A  zero  indicate  a  missile  launched  at  time 
interval  u  is  not  airborne  during  time  period  t.  An 
additional  five  seconds  is  added  to  each  time  of  flight 
to  allow  for  operator  assessment  of  the  engagement. 

3 .  Decision  Variable 


X^:  A  binary  variable,  a  one  indicates  a  missile  should  be 
launched  during  time  interval  t,  a  zero  indicates  no 
missile  is  to  be  fired  dxuring  time  interval  t.  Although 
is  a  binary  variable,  in  practice  it  only  needs  to  be 
limited  to  values  less  than  or  equal  to  one  to  provide 
binary  results. 


4 .  Foxnulatlon 


Minimize  ^ 
e-i 


(Obj) 


Subject  to: 

VI 

(1) 

2}  ^  2,  foz  all  u 

C«1 

(2) 

i  1 ,  a  tiMv 

(3) 

Xf  1,  for  all  t 

(4) 

In  these  equations,  constraint  (1)  limits  the  number  of 
missiles  fired  to  a  maximum  of  4  per  engagement,  constraint 
(2)  limits  the  number  of  missiles  airborne  at  any  time  to  two, 
constraint  (3)  limits  the  rate  of  fire  to  one  loissile  every 
two  seconds,  and  constraint  (4)  limits  the  number  of  missiles 
fired  during  any  given  time  period  to  one. 

B.  BXBDLTS 

To  detejnnine  the  effect  of  using  the  optimized  firing 
doctrine,  the  linear  program  was  in^lemented  using  the  General 
Algebraic  Modeling  System  (GAMS) .  To  avoid  using  classified 
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data,  a  range  of  ten  nautical  miles  (NM)  and  a  constant 
velocity  of  600  nautical  miles  per  hour  (kts)  was  assumed  for 
the  Seasparrow.  The  Pk  versus  Range  curve  found  in  figures  1 
and  2  was  utilized,  and  the  Exocet  was  assumed  to  have  a 
constant  speed  of  600  kts-  The  optimized  firing  doctrine 
returned  a  value  for  Pk(4)j  of  0.967,  A  shoot -shoot -look- 
shoot-shoot  doctrine  with  the  first  shots  at  near  maximum 
range  returns  a  value  of  0.927,  a  difference  of  0.04.  An  even 
bigger  difference  is  seen  when  the  probability  of  destroying 
the  target  with  the  first  two  missiles  is  examined.  With  a 
shoot  when  the  target  is  within  range  doctrine,  the 
probability  of  destroying  the  target  is  0.584.  The  optimized 
doctrine  increases  the  probability  of  success  to  0.752.  This 
is  especially  significant  given  the  limited  number  of  missiles 
available  in  the  launcher  and  the  extended  time  required  to 
reload  the  launcher.  In  all  cases,  the  increased  Pk  from  the 
optimized  doctrine  results  from  a  shorter  intercept  range, 
with  a  corresponding  increase  in  single  round  Pk. 

Although  these  numbers  are  only  for  illustrative  purposes, 
they  show  the  significant  increase  available  from  an  optimized 
firing  doctrine  when  compared  to  a  doctrine  that  assumes  a 
constant  Pk-  The  optimization  prograuti  ran  on  a  80486/ 3 3MH2 
con^uter  in  just  tinder  1.5  seconds,  and  can  be  easily  modified 
to  accommodate  almost  any  firing  doctrine  by  sinply  changing 
constants.  The  GAMS  model  is  provided  in  Appendix  B.  Before 
using  with  actual  data,  the  missile  in-flight  confotinding 
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matrix,  C(u,t)  ,  needs  to  be  updated  with  actual  time  of  flight 
data. 
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VIZ.  CONCLUSIONS  AMD  RgCOMMKNPATIOMS 


Dr.  Woods'  baseline  system  effectiveness  model  for  the 
NSSMS  provides  the  framework  to  assess  both  the  availability 
and  capability  of  the  system.  This  thesis  provides  em. 
enhancement  to  the  baseline  model  by  providing  a  method  to 
determine  the  probability  of  kill  for  the  NATO  Seasparrow 
missile  as  a  function  of  range  and  system  state.  The  Pk  versus 
range  curve  can  be  developed  using  existing  missile 
simulations  and,  when  used  in  conjunction  with  the  available 
inissile  reliability  data,  can  provide  a  way  to  estimate  the 
capability  of  the  NSSMS  to  defeat  the  designated  threat. 

Beyond  enhancing  the  system  effectiveness  model,  a  Pk 
versus  range  curve  can  provide  insight  into  the  proper 
tactical  enployment  of  the  system.  Through  the  application  of 
linear  programming  techniques,  it  is  possible  to  significantly 
increase  the  probability  of  success  by  changing  the  firing 
doctrine  to  irtaximize  Pk.  Additionally,  the  development  of  Pk 
versus  range  curves  for  degraded  states  of  operation  cam 
provide  the  basis  for  firing  strategies  to  minimize  the 
effects  of  the  degradations. 

Further  enhancements  to  the  baseline  system  effectiveness 
model  should  concentrate  on  incorporating  the  effects  of 
search  radars  and  acquisition  times  into  the  model,  as 
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outlined  in  chapter  V.  The  baseline  model  could  easily  be 
expanded  to  include  NSSHS,  search  radars,  and  close  in  weapons 
systems  (CIWS) ,  providing  an  assessment  o£  the  capability  of 
a  ship  to  defend  itself  against  a  chosen  threat. 

Modem  weapon  systems  with  several  degraded  modes  of 
operation  make  the  traditional  methods  of  tracking  systems 
availability  less  suitable.  The  baseline  systems  effectiveness 
model  for  the  NSSMS,  with  the  enhemcements  provided  by  this 
paper,  provides  the  technical  community  with  a  better  method 
of  measuring  the  utility  of  the  NATO  Seasparrow. 
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APPENDIX  A.  SniDLATZCW  PLAN 


The  following  simulation  plan  is  proposed  for  execution  by  the 
Naval  Air  Warfare  Center,  China  Lake,  California,  using  the  six 
degree  of  freedom,  digital  Sparrow  simulation  and  associated 
endgeune  lethality  simulations.  The  procedure  for  requesting 
utilization  of  simulations  in  the  NSPO  performance  assessment 
models  network  to  conduct  studies  is  outlined  in  Reference  1. 


Common  to  all  cases. 

Missile;  RIM-7P 
Target;  Exocet  SSM 

Target  Profile;  Exocet  profile  with  firing  ship  as  the  intended 
target.  The  Exocet  acquires  the  firing  ship  at  laaximum  range 
and  commences  sea  skimming  at  the  earliest  opportunity. 

Simulation  Points; 

Intercept  Remae  (%  of  maximum)  Number  of  runs 


5.0 

50 

7.5 

50 

10.0 

50 

15.0 

50 

20.0 

50 

30.0 

50 

40.0 

50 

55.0 

50 

70.0 

50 

85.0 

50 

90.0 

50 

95.0 

50 

100.0 

50 

All  range  at  launch  missile  command  degradations  are 
created  by  adjusting  the  actual  launch  range,  not  the  range  at 
launch  commands.  The  range  at  launch  command  should  always  reflect 
the  range  required  for  intercept  at  the  desired  range. 
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Case  1;  NSSMS  system  states  Sj :  All  up,  S^:  FOC/RSC  Degraded,  and 
S^:  Conputer  Complex  Casualty  Mode. 

Missile  Commands:  Nominal 

Launcher  Super-Elevation:  In  accordance  with  algorithm 
Transmitter  Power:  Nominal 


Case  2:  NSSMS  system  state  S2:  Transmitter  power  degraded. 
Missile  Commands;  Nominal 

Launcher  Super-Elevation:  In  accordance  with  algorithm 
Transmitter  Power:  Degraded 


Case  3;  NSSMS  system  state  S3:  SDP  range  signal  degraded. 

Missile  Commands;  Rzuige  at  launch  deviation  i6%-15%  from 
actual,  all  others  nominal 

Launcher  Super-Elevation:  In  accordance  with  algorithm 
Transmitter  Power:  Nominal 


Case  4:  NSSMS  system  state  S^:  SDP  Range  deviation  too  large. 

Missile  Commands:  Ramge  at  launch  deviation  greater  than 
+15%  from  actual,  all  others  nominal 
Launcher  Super-Elevation:  In  accordance  with  algorithm 
Transmitter  Power:  Nominal 


Case  5:  NSSMS  system  state  Sg:  Slaved  Illuminator/Low  Light 
Level  Television  mode. 

Missile  Commands:  Default 
Launcher  Super-Elevation;  Default 
Transmitter  Power;  Nominal 


Note:  Values  for  missile  command  degradations  presented  here 
and  in  Reference  1,  Appendix  F  are  tentative,  and  will  require 
modification  following  the  NATO  Seasparrow  technical  community 
review  of  Reference  1. 
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APPENDIX  B.  GAMS  SOORCE  OQOE  FCR  (XTIMIZB>  DOCIRINE 


$TITLE  Thesis  Formulatitai 

$STITLE  Optimized  Firing  Doctrine 

* - gams  and  dollar  control  options 


$OTFUPPER  OFFSYMLIST  OFFSXMXREF 

options 

LIMCOL  =  0  ,  LIMROW  =  0  ,  SOLPRINT  =  OFF  ,  DECIMALS  =  2 

RESLIM  =  100  ,  ITERLIM  =  10000  ,  OPTCR  =  0.00  ,  SEED  =  3141; 

* - DEFINITIOTS  AND  DATA - 

SCALAR  TM  total  missiles  available  /4/ 

AM  missiles  airborne  at  one  time  /2/; 

SETS 

t  lamch  time  interval 
tl  subset  of  t 

t2  subset  of  t 

t3  subset  of  t 

t4  subset  of  t 

tS  subset  of  t 

PARAMETERS 


PKl(tl)  negative 

In  of  1-: 

71 

0.2881 

2 

0.3173 

3 

0.3445 

4 

0.3697 

5 

0.3927 

6 

0.4137 

7 

0.4326 

8 

0.4495 

9 

0.4646 

10 

0.4779 

11 

0.4897 

12 

0.5001 

13 

0.5093 

14 

0.5173 

15 

0.5245 

16 

0.5310 

17 

0.5368 

18 

0.5422 

19 

0.5472 

20 

0.5520 

21 

0.5566 

/1*115/ 

71*25/ 

726*50/ 

751*75/ 

776*100/ 

7101*1157; 
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22 

0.5612 

23 

0.5657 

24 

0.5702 

25 

0.5748/ 

PK2(t2)  negative 

In  of  1-Pk 

/26 

0.5794 

27 

0.5842 

28 

0.5890 

29 

0.5939 

30 

0.5989 

31 

0.6040 

32 

0.6091 

33 

0.6142 

34 

0.6194 

35 

0.6245 

36 

0.6296 

37 

0.6347 

38 

0.6397 

39 

0.6446 

40 

0.6494 

41 

0.6541 

42 

0.6588 

43 

0.6633 

44 

0.6677 

45 

0.6720 

46 

0.6763 

47 

0.6804 

48 

0.6846 

49 

0.6887 

50 

0.6928/ 

PK3(t3)  negative 

In  of  1-Pk 

/51 

0.6969 

52 

0.7011 

53 

0.7054 

54 

0.7097 

55 

0.7142 

56 

0.7188 

57 

0.7236 

58 

0.7285 

59 

0.7336 

60 

0.7389 

61 

0.7445 

62 

0.7502 

63 

0.7560 

64 

0.7621 

65 

0.7683 

66 

0.7747 

67 

0.7811 

68 

0.7876 

69 

0.7942 

70 

0.8007 

35 


71 

0.8072 

72 

0.8136 

73 

0.8199 

74 

0.8259 

75 

0.8317/ 

PK4(t4)  negative  In  of  1-Pk 

/76 

0.8373 

77 

0.8425 

78 

0.8474 

79 

0.8519 

80 

0.8560 

81 

0.8597 

82 

0.8631 

83 

0.8661 

84 

0.8688 

85 

0.8711 

86 

0.8733 

87 

0.8753 

88 

0.8771 

89 

0.8789 

90 

0.8808 

91 

0.8828 

92 

0.8850 

93 

0.8874 

94 

0.8901 

95 

0.8930 

96 

0.8961 

97 

0.8993 

98 

0.9023 

99 

0.9050 

100 

0.9070/ 

PK5(t5)  negative 

In  of  1-Pk 

/lOl 

0.9076 

102 

0.9063 

103 

0.9021 

104 

0.8941 

105 

0.8812 

106 

0.8618 

107 

0.8348 

108 

0.7985 

109 

0.7518 

110 

0.6935 

111 

0.6229 

112 

0.5396 

113 

0.4439 

114 

0.3364 

115 

0.2181/; 

TABLE 

Cl(t,tl)  nvissile  in  air  confounding  natrix 

1  2  3  4  5  6  7  8  9  10  11  12  13  14  15  16  17  18  19  20  21  22  23  24  25 
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1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

2 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

3 

0 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

4 

0 

0 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

5 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

6 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

7 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

8 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

9 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

10 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

11 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

12 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

13 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

14 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

15 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

16 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

1 

1 

1 
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0 

0 

0 
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0 
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0 
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0 
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54  00000000 

55  00000000 
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59  00000000 

60  00000000 
61  00000000 
62  00000000 

63  00000000 

64  00000000 

65  00000000 

66  00000000 

67  0  0  0  0  0  0  0  0 

68  0  0  0  0  0  0  0  0 

69  00000000 

70  00000000 
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86  00000000 
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91  00000000 
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95  0  0  0  0  0  0  0  0 

96  00000000 
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99  00000000 

100  00000000 
101  00000000 
102  00000000 

103  00000000 

104  00000000 

105  00000000 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

113 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0; 

TABLE 

C2(t,t2)  missile  in  air  confounding  matrix 

26  27  28  29  30  31  32  33  34  35  36  37  38  39  40  41  42  43  44  45  46  47  48  49  50 
1  1111111111111111111111111 

2  1111111111111111111111111 

3  1111111111111111111111111 

4  1111111111111111111111111 

5  1111111111111111111111111 

6  1111111111111111111111111 

7  1111111111111111111111111 

8  1111111111111111111111111 

9  1111111111111111111111111 

10  1111111111111111111111111 

11  1111111111111111111111111 

12  1111111111111111111111111 

13  1111111111111111111111111 

14  1111111111111111111111111 

15  1111111111111111111111111 

16  1111111111111111111111111 

17  1111111111111111111111111 

18  1111111111111111111111111 

19  1111111111111111111111111 

20  1111111111111111111111111 

21  1111111111111111111111111 

22  1111111111111111111111111 

23  1111111111111111111111111 

24  1111111111111111111111111 

25  1111111111111111111111111 

26  1111111111111111111111111 

27  0111111111111111111111111 

28  0011111111111111111111111 

29  0001111111111111111111111 

30  0000111111111111111111111 

31  0000011111111111111111111 

32  0000001111111111111111111 

33  0000000111111111111111111 

34  0000000011111111111111111 

35  0000000001111111111111111 

36  0000000000111111111111111 

37  0000000000011111111111111 

38  0000000000001111111111111 

39  0000000000000111111111111 


39 


40 


93  0000000000000000000000000 

94  0000000000000000000000000 

95  0000000000000000000000000 

96  0000000000000000000000000 

97  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0 

98  0000000000000000000000000 

99  0000000000000000000000000 

100  0000000000000000000000000 

101  0000000000000000000000000 

102  0000000000000000000000000 

103  0000000000000000000000000 

104  0000000000000000000000000 

105  0000000000000000000000000 

106  0000000000000000000000000 

107  0000000000000000000000000 

108  0000000000000000000000000 

109  0000000000000000000000000 

110  0000000000000000000000000 

111  0000000000000000000000000 

112  0000000000000000000000000 

113  0000000000000000000000000 

114  0000000000000000000000000 

115  000000000000000000000000  0; 


TABLE 

C3(t,t3)  nvLssile  in  air  confounding  matrix 

51  52  53  54  55  56  57  58  59  60  61  62  63  64  65  66  67  68  69  70  71  72  73  74  75 
1  1111111111111110000000000 

2  1111111111111111000000000 

3  1111111111111111000000000 

4  1111111111111111100000000 

5  1111111111111111100000000 

6  1111111111111111110000000 

7  1111111111111111110000000 

8  1111111111111111111000000 

9  1111111111111111111000000 

10  1111111111111111111100000 

11  1111111111111111111100000 

12  1111111111111111111110000 

13  1111111111111111111110000 

14  1111111111111111111111000 

15  1111111111111111111111000 

16  1111111111111111111111100 

17  1111111111111111111111100 

18  1111111111111111111111110 

19  1111111111111111111111110 

20  1111111111111111111111111 

21  1111111111111111111111111 

22  1111111111111111111111111 

23  1111111111111111111111111 

24  1111111111111111111111111 

25  1111111111111111111111111 


41 


26  11111111111 

27  11111111111 

28  11111111111 

29  11111111111 

30  11111111111 

31  11111111111 

32  11111111111 

33  11111111111 

34  11111111111 

35  11111111111 

36  11111111111 

37  11111111111 

38  11111111111 

39  11111111111 

40  11111111111 

41  11111111111 

42  11111111111 

43  11111111111 

44  11111111111 

45  11111111111 

46  11111111111 

47  11111111111 

48  11111111111 

49  11111111111 

50  11111111111 

51  11111111111 

52  01111111111 

53  00111111111 

54  00011111111 

55  00001111111 

56  00000111111 

57  00000011111 

58  00000001111 

59  00000000111 

60  0  0  0  0  0  0  0  0  011 

61  00000000001 

62  00000000000 

63  00000000000 

64  00000000000 

65  00000000000 

66  00000000000 

67  0  0  0  0  0  0  0  0  0  0  0 

68  00000000000 

69  00000000000 

70  00000000000 

71  00000000000 

72  00000000000 

73  00000000000 

74  00000000000 

75  00000000000 

76  00000000000 

77  00000000000 

78  00000000000 


11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

11111111111111 

01111111111111 

00111111111111 

00011111111111 

00001111111111 

00000111111111 

00000011111111 

00000001111111 

00000000111111 

00000000011111 

00000000001111 

00000000000111 

00000000000011 

00000000000001 

00000000000000 

00000000000000 

00000000000000 


42 


79  0000000000000000000000000 

80  0000000000000000000000000 

81  0000000000000000000000000 

82  0000000000000000000000000 

83  0000000000000000000000000 

84  0000000000000000000000000 

85  0000000000000000000000000 

86  0000000000000000000000000 

87  0000000000000000000000000 

88  0000000000000000000000000 

89  0000000000000000000000000 

90  0000000000000000000000000 

91  0000000000000000000000000 

92  0000000000000000000000000 

93  0000000000000000000000000 

94  0000000000000000000000000 

95  0000000000000000000000000 

96  0000000000000000000000000 

97  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0  0 

98  0000000000000000000000000 

99  0000000000000000000000000 

100  0000000000000000000000000 

101  0000000000000000000000000 

102  0000000000000000000000000 

103  0000000000000000000000000 

104  0000000000000000000000000 

105  0000000000000000000000000 

106  0000000000000000000000000 

107  0000000000000000000000000 

108  0000000000000000000000000 

109  0000000000000000000000000 

110  0000000000000000000000000 

111  0000000000000000000000000 

112  0000000000000000000000000 

113  0000000000000000000000000 

114  0000000000000000000000000 

115  000000000000000000000000  0; 

TABLE 

C4(t,t4)  missile  in  air  confounding  matrix 

76  77  78  79  80  81  82  83  84  85  86  87  88  89  90  91  92  93  94  95  96  97  98  99  100 

1  000000000000000000000000  0 

2  000000000000000000000000  0 

3  000000000000000000000000  0 

4  000000000000000000000000  0 

5  000000000000000000000000  0 

6  000000000  0  00000000000000  0 

7  000000000000000000000000  0 

8  000000000000000000000000  0 

9  000000000000000000000000  0 

10  000000000000000000000000  0 

11  000000000000000000000000  0 

12  000000000000000000000000  0 
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66  11111111111 

67  11111111111 

68  11111111111 

69  11111111111 

70  11111111111 

71  11111111111 

72  11111111111 

73  11111111111 

74  11111111111 

75  11111111111 

76  11111111111 

77  01111111111 

78  00111111111 

79  00011111111 

80  00001111111 

81  00000111111 

82  00000011111 

83  00000001111 

84  00000000111 

85  00000000011 

86  00000000001 

87  00000000000 

88  000000000000 

89  000000000000 

90  0  0  0  0  0  0  0  0  0  0  0  0 

91  000000000000 

92  000000000000 

93  000000000000 

94  000000000000 

95  000000000000 

96  000000000000 

97  000000000000 

98  000000000000 

99  000000000000 

100  000000000000 

101  000000000000 

102  000000000000 

103  000000000000 

104  000000000000 

105  000000000000 

106  000000000000 

107  000000000000 

108  00000  000000 

109  000000000000 

110  000000000000 

111  000000000000 

112  000000000000 

113  000000000000 

114  oooooooooooc 

115  000000000000 

TABLE 

C5(t,t5)  missile  in  air  confounding  matrix 


1111111111111  0 

1111111111111  0 

1111111111111  1 

1111111111111  1 

1111111111111  1 

1111111111111  1 

1111111111111  1 

1111111111111  1 

1111111111111  1 

1111111111111  1 

1111111111111  1 

1111111111111  1 

1111111111111  1 

1111111111111  1 

1111111111111  1 

1111111111111  1 

1111111111111  1 

1111111111111  1 

1111111111111  1 

111111111111  1 

111111111111  1 

111111111111  1 

111111111111  1 

011111111111  1 

001111111111  1 

000111111111  1 

000011111111  1 

000001111111  1 

000000111111  1 

000000011111  1 

000000001111  1 

O'^OOOOOOOlll  1 

0  0000000011  1 

0  )0000000001  1 

000000000000  1 

000000000000  0 

000000000000  0 

000000000000  0 

000000000000  0 

000000000000  0 

000000000000  0 

000000000000  0 

000000000000  0 

000000000000  0 
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105  00001111 

106  00000111 

107  00000011 

108  00000001 
109  00000000 

no  00000000 
111  00000000 
112  00000000 

113  00000000 

114  00000000 

115  00000000 


TABLE 

Rl(t,tl)  rate  of  fire  limitations 
1  2  3  4  5  6  7  8  9  10  11  12  13 

1  110000000  0  0  0  0 

2  011000000  0  0  0  0 

3  001100000  0  0  0  0 

4  000110000  0  0  0  0 

5  000011000  0  0  0  0 

6  000001100  0  0  0  0 

7  000000110  0  0  0  0 

8  000000011  0  0  0  0 

9  000000001  1  0  0  0 

10  000000000  1  1  0  0 

11  000000000  0  1  1  0 

12  0000000000011 

13  0000000000001 

14  0000000000000 

15  000000000  0  0  0  0 

16  000000000  0  0  0  0 

17  000000000  0  0  0  0 

18  0000000000000 

19  000000000  0  0  0  0 

20  000000000  0  0  0  0 

21  000000000  0  0  0  0 

22  0000000000000 

23  000000000  0  0  0  0 

24  000000000  0  0  0  0 

25  0000000000000 

26  000000000  0  0  0  0 

27  000000000  0  0  0  0 

28  0000000000000 

29  000000000  0  0  0  0 

30  000000000  0  0  0  0 

31  000000000  0  0  0  0 

32  000000000  0  0  0  0 

33  0000000000000 

34  0000000000000 

35  0000000000000 

36  0000000000000 

37  000000000  0  0  0  0 


1111111 
1111111 
1111111 
1111111 
1111111 
0  111111 

0  0  11111 

0  0  0  1  1  1  1 

0  0  0  0  1  1  1 

0  0  0  0  0  1  1 

0  0  0  0  0  0  1; 


14  15  16  17  18  19  20  21  22  23  24  25 
000000000000 
000000000000 
000000000000 
000000000000 
000000000000 
000000000000 
000000000000 
000000000000 
000000000000 
000000000000 
000000000000 
000000000000 
100000000000 
110000000000 
011000000000 
001100000000 
000110000000 
000011000000 
000001100000 
000000110000 
000000011000 
000000001100 
000000000110 
000000000011 
000000000001 
000000000000 
000000000000 
000000000000 
000000000000 
000000000000 
000000000000 
000000000000 
000000000000 
000000000000 
000000000000 
000000000000 
000000000000 
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38  000000000  0  0  0  0 

39  000000000  0  0  0  0 

40  000000000  0  0  0  0 

41  000000000  0  0  0  0 

42  000000000  0  0  0  0 

43  000000000  0  0  0  0 

44  000000000  0  0  0  0 

45  0000000000000 

46  0000000000000 

47  000000000  0  0  0  0 

48  000000000  0  0  0  0 

49  0000000000000 

50  000000000  0  0  0  0 

51  0000000000000 

52  000000000  0  0  0  0 

53  000000000  0  0  0  0 

54  0000000000000 

55  000000000  0  0  0  0 

56  000000000  0  0  0  0 

57  0  0  0  0  0  0  0  0  0  0  0  0  0 

58  000000000  0  0  0  0 

59  0000000000000 

60  000000000  0  0  0  0 

61  000000000  0  0  0  0 

62  000000000  0  0  0  0 

63  000000000  0  0  0  0 

64  000000000  0  0  0  0 

65  000000000  0  0  0  0 

66  0000000000000 

67  000000000  0  0  0  0 

68  000000000  0  0  0  0 

69  0000000000000 

70  0000000000000 

71  000000000  0  0  0  0 

72  000000000  0  0  0  0 

73  000000000  0  0  0  0 

74  000000000  0  0  0  0 

75  000000000  0  0  0  0 

76  000000000  0  0  0  0 

77  0000000000000 

78  000000000  0  0  0  0 

79  000000000  0  0  0  0 

80  000000000  0  0  0  0 

81  0000000000000 
82  0000000000000 

83  0000000000000 

84  000000000  0  0  0  0 

85  0000000000000 

86  000000000  0  0  0  0 

87  000000000  0  0  0  0 

88  0000000000000 

89  000000000  0  0  0  0 

90  0  0  0  0  0  0  0  0  0  0  0  0  0 


000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 
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000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 

000000000000 


49 


91 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

92 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

93 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

94 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

95 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

96 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

97 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

98 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

99 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

100 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

101 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

’.02 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

103 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

104 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

105 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

106 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

107 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

108 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

109 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

110 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

111 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

112 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

113 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

114 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

115 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0; 

TABLE 

R2{t,t2)  rate  of  fire  limitations 

26  27  28  29  30  31  32  33  34  35  36  37  38  39  40  41  42  43  44  45  46  47  48  49  50 
1  0000000000000000000000000 

2  0000000000000000000000000 

3  0000000000000000000000000 

4  0000000000000000000000000 

5  0000000000000000000000000 

6  0000000000000000000000000 

7  0000000000000000000000000 

8  0000000000000000000000000 

9  0000000000000000000000000 

10  0000000000000000000000000 

11  0000000000000000000000000 

12  0000000000000000000000000 

13  0000000000000000000000000 

14  0000000000000000000000000 

15  0000000000000000000000000 

16  0000000000000000000000000 

17  0000000000000000000000000 

18  0000000000000000000000000 

19  0000000000000000000000000 

20  0000000000000000000000000 

21  0000000000000000000000000 

22  0000000000000000000000000 
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23  00000000 

24  00000000 

25  10000000 

26  1  1000000 

27  01100000 

28  00110000 

29  00011000 

30  00001100 

31  00000110 

32  00000011 

33  00000001 

34  00000000 

35  00000000 

36  00000000 

37  00000000 

38  00000000 

39  00000000 

40  00000000 

41  00000000 

42  00000000 

43  00000000 

44  00000000 

45  00000000 

46  00000000 

47  00000000 

48  00000000 

49  00000000 

50  00000000 

51  00000000 

52  00000000 

53  00000000 

54  00000000 

55  00000000 

56  00000000 

57  00000000 

58  00000000 

59  00000000 

60  00000000 
61  00000000 
62  00000000 

63  00000000 

64  00000000 

65  00000000 

66  00000000 

67  00000000 

68  00000000 

69  00000000 

70  00000000 

71  00000000 

72  00000000 

73  00000000 

74  00000000 

75  00000000 


00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

10000000 

11000000 

01100000 

00110000 

00011000 

00001100 

00000110 

00000011 

00000001 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 

00000000 


000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

100000000 

110000000 

011000000 

001100000 

000110000 

000011000 

000001100 

000000110 

000000011 

000000001 

000300000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 
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76  000000000000 

77  000000000000 

78  000000000000 

79  000000000000 

80  000000000000 
81  000000000000 
82  000000000000 

83  000000000000 

84  000000000000 

85  000000000000 

86  000000000000 

87  000000000000 

88  000000000000 

89  000000000000 

90  000000000000 

91  000000000000 

92  000000000000 

93  000000000000 

94  000000000000 

95  000000000000 

96  000000000000 

97  000000000000 

98  000000000000 

99  000000000000 

100  000000000000 
101  000000000000 
102  000000000000 

103  000000000000 

104  000000000000 

105  000000000000 

106  000000000000 

107  000000000000 

108  000000000000 

109  000000000000 

110  000000000000 
111  000000000000 
112  000000000000 

113  000000000000 

114  000000000000 

115  000000000000 


0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
0000000000000 
000000000000  0; 


TABLE 

R3(t,t3)  rate  of  fire  limitations 

51  52  53  54  55  56  57  58  59  60  61  62  63 
1  0000000000000 

2  0000000000000 

3  0000000000000 

4  0000000000000 

5  0000000000000 

6  0000000000000 

7  0000000000000 

8  0000000000000 

9  0000000000000 


64  65  66  67 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 


68  69  70  71 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 


72  73  74  75 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
0  0  0  0 
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10 
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0 

0 

0 
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0 

0 

0 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

11 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

12 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

13 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

14 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

15 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

16 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

17 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

18 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

19 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

20 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

21 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

22 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

23 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

24 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

25 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

26 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

27 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

28 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

29 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

30 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

31 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

32 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

33 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

34 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

35 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

36 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

37 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

38 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

39 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

40 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

41 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

42 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

43 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

44 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

45 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

46 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

47 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

48 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

49 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

50 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

51 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

52 

0 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

53 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

54 

0 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

55 

0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

56 

0 

0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

57 

0 

0 

0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

58 

0 

0 

0 

0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

59 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

60 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

61 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

62 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

53 


63  0  0  0  0  0  0 

64  0  0  0  0  0  0 

65  0  0  0  0  0  0 

66  0  0  0  0  0  0 

67  0  0  0  0  0  0 

68  0  0  0  0  0  0 

69  0  0  0  0  0  0 

70  0  0  0  0  0  0 

71  0  0  0  0  0  0 

72  0  0  0  0  0  0 

73  0  0  0  0  0  0 

74  0  0  0  0  0  0 

75  0  0  0  0  0  0 

76  0  0  0  0  0  0 

77  0  0  0  0  0  0 

78  0  0  0  0  0  0 

79  0  0  0  0  0  0 

80  0  0  0  0  0  0 

81  0  0  0  0  0  0 

82  0  0  0  0  0  0 

83  0  0  0  0  0  0 

84  0  0  0  0  0  0 

85  0  0  0  0  0  0 

86  0  0  0  0  0  0 

87  0  0  0  0  0  0 

88  0  0  0  0  0  0 

89  0  0  0  0  0  0 

90  0  0  0  0  0  0 

91  0  0  0  0  0  0 

92  0  0  0  0  0  0 

93  0  0  0  0  0  0 

94  0  0  0  0  0  0 

95  0  0  0  0  0  0 

96  0  0  0  0  0  0 

97  0  0  0  0  0  0 

98  0  0  0  0  0  0 

99  0  0  0  0  0  0 

100  0  0  0  0  0  0 

101  0  0  0  0  0  0 

102  0  0  0  0  0  0 

103  0  0  0  0  0  0 

104  0  0  0  0  0  0 

105  0  0  0  0  0  0 

106  0  0  0  0  0  0 

107  0  0  0  0  0  0 

108  0  0  0  0  0  0 

109  0  0  0  0  0  0 

110  0  0  0  0  0  0 

111  0  0  0  0  0  0 

112  0  0  0  0  0  0 

113  0  0  0  0  0  0 

114  0  0  0  0  0  0 

115  0  0  0  0  0  0 


000000110 

000000011 

000000001 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

uOuOOOOGf^ 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 


0000000000 

0000000000 

1000000000 

1100000000 

0110000000 

0011000000 

0001100000 

0000110000 

0000011000 

0000001100 

0000000110 

0000000011 

0000000001 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000 

ooooocoooo 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000 

0000000000; 


54 


TABLE 

R4(t,t4)  rate  of  fire  limitations 

76  77  78  79  80  81  82  83  84  85  86 
1  00000000000 

2  00000000000 

3  00000000000 

4  00000000000 

5  00000000000 

6  00000000000 

7  00000000000 

8  00000000000 

9  00000000000 

10  00000000000 

11  00000000000 

12  00000000000 

13  00000000000 

14  00000000000 

15  00000000000 

16  00000000000 

17  00000000000 

18  00000000000 

19  00000003000 

20  00000000000 

21  00000000000 

22  00000000000 

23  00000000000 

24  00000000000 

25  00000000000 

26  00000000000 

27  00000000000 

28  00000000000 

29  00000000000 

30  00000000000 

31  00000000000 

32  00000000000 

33  00000000000 

34  00000000000 

35  00000000000 

36  00000000000 

37  00000000000 

38  00000000000 

39  00000000000 

40  00000000000 

41  00000000000 

42  00000000000 

43  00000000000 

44  00000000000 

45  00000000000 

46  00000000000 

47  00000000000 

48  00000000000 

49  00000000000 


87  88  89  90  91  92  93  94  95  96  97  98  99  100 
0000000000000  0 
0000000000000  0 
0000000000000  0 
0000000000000  0 
0000000000000  0 
0000000000000  0 
0000000000000  0 
0000000000000  0 
0000000000000  0 
0000000000000  0 
0000000000000  0 
0000000000000  0 
0000000000000  0 
0000000000000  0 
0000000000000  0 
0000000000000  0 
0000000000000  0 
0000000000000  0 
0000000000000  0 
0000000000000  0 
0000000000000  0 
0000000000000  0 
OOOOOOOOOOOOO  0 
0000000000000  0 
OOOOOOOOOOOOO  0 
OOOOOOOOOOOOO  0 
OOOOOOOOOOOOO  0 
OOOOOOOOOOOOO  0 
OOOOOOOOOOOOO  0 
OOOOOOOOOOOOO  0 
OOOOOOOOOOOOO  0 
OOOOOOOOOOOOO  0 
OOOOOOOOOOOOO  0 
OOOOOOOOOOOOO  0 
OOOOOOOOOOOOO  0 
OOOOOOOOOOOOO  0 
OOOOOOOOOOOOO  0 
OOOOOOOOOOOOO  0 
OOOOOOOOOOOOO  0 
OOOOOOOOOOOOO  0 
OOOOOOOOOOOOO  0 
OOOOOOOOOOOOO  0 
OOOOOOOOOOOOO  0 
OOOOOOOOOOOOO  0 
OOOOOOOOOOOOO  0 
OOOOOOOOOOOOO  0 
OOOOOOOOOOOOO  0 
OOOOOOOOOOOOO  0 
OOOOOOOOOOOOO  0 


55 


50  0  0  0  0  0  0 

51  0  0  0  0  0  0 

52  0  0  0  0  0  0 

53  0  0  0  0  0  0 

54  0  0  0  0  0  0 

55  0  0  0  0  0  0 

56  0  0  0  0  0  0 

57  0  0  0  0  0  0 

58  0  0  0  0  0  0 

59  0  0  0  0  0  0 

60  0  0  0  0  0  0 

61  0  0  0  0  0  0 

62  0  0  0  0  0  0 

63  0  0  0  0  0  0 

64  0  0  0  0  0  0 

65  0  0  0  0  0  0 

66  0  0  0  0  0  0 

67  0  0  0  0  0  0 

68  0  0  0  0  0  0 

69  0  0  0  0  0  0 

70  0  0  0  0  0  0 

71  0  0  0  0  0  0 

72  0  0  0  0  0  0 

73  0  0  0  0  0  0 

74  0  0  0  0  0  0 

75  1  0  0  0  0  0 

76  1  1  0  0  0  0 

77  0  1  1  0  0  0 

78  0  0  1  1  0  0 

79  0  0  0  1  1  0 

80  0  0  0  0  1  1 

81  0  0  0  0  0  1 

82  0  0  0  0  0  0 

83  0  0  0  0  0  0 

84  0  0  0  0  0  0 

85  0  0  0  0  0  0 

86  0  0  0  0  0  0 

87  0  0  0  0  0  0 

88  0  0  0  0  0  0 

89  0  0  0  0  0  0 

90  0  0  0  0  0  0 

91  0  0  0  0  0  0 

92  0  0  0  0  0  0 

93  0  0  0  0  0  0 

94  0  0  0  0  0  0 

95  0  0  0  0  0  0 

96  0  0  0  0  0  0 

97  0  0  0  0  0  0 

98  0  0  0  0  0  0 

99  0  0  0  0  0  0 

100  0  0  0  0  0  0 

101  0  0  0  0  0  0 

102  0  0  0  0  0  0 


000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

100000000 

110000000 

011000000 

001100000 

000110000 

000011000 

000001100 

000000110 

000000011 

000000001 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 


000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

000000000 

100000000 

110000000 

011000000 

001100000 

000110000 

000011000 

000001100 

000000110 

000000011 

000000001 

000000000 

000000000 

000000000 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1 

0 

0 


56 


103 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

104 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

105 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

106 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

107 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

108 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

109 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

110 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

111 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

112 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

113 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

114 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

115 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0; 

TABLE 

R5(t,t5) 
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■MODEL 


BINARY  VARIABLES 

X(t)  indicates  if  firing  in  time  period 
Xl(tl) 

X2(t2) 

X3(t3) 

X4(t4) 

X5(t5); 


VARIABLE 

TPK  total  pk; 


EQUATICWS 

OBJ  total  pk 

NMIS  number  of  missiles  available 

AMIS(t)  number  of  missiles  airborne  at  once 

ROF(t)  rate  of  fire  limitations; 


*  MAXIMIZE 
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OBJ..  TPK  =E=  SUM(tl,  Xl(tl)*PKl(tl))+SUM(t2,  X2(t2)*PK2(t2) ) 
+SUM(t3,  X3(t3)*PK3(t3))+SUM(t4,  X4(t4)*PK4(t4) ) 
+SUM(t5,  X5(t:')*PK5(t5)); 

*  subject  to 

NMIS..  SUM(tl,  Xl(tl))+SUM(t2,  X2(t2))+SUM(t3,  X3(t3)) 
+SUM(t4,  X4(t4))+SUM(t5,  X5(t5))  =L=  IM; 

AMIS(t)..  SUM(tl,  Cl(t,tl)*Xl(fcl))+SUM(t2,  C2(t ,t2)*X2(t2) ) 
+SUM(t3,  C3(t,t3)*X3(t3))+SUM(t4  ,C4(t ,t4)*X4(t4) ) 
+SUMU5  ,C5(t.t5)*X5(t5))  =L=  Wd; 

ROF(t)..  SUM(tl,  Rl(t,tl)*Xl(tl))+SUM(t2,  R2(t,t2)*X2(t2)) 
+SUM(t3,  R3(t,t3)*X3(t3))+SUM(t4,  R4(t,t4)*X4(t4)) 
+SUM(t5,  R5(t,t5)*X5(t5))  =L=  1; 

MODEL  MIX  /ALL/ 

SOLVE  MIX  USING  MIP  MAXIMIZING  TPK; 

*  - REPORTS - 

*  print  the  optimal  objective  value  and  solution. 

DISPLAY  TPK.L,  XI. L,  X2.L,  X3.L,  X4.L,  X5.L; 
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